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The product of the branching ratios of the standard model (SM) Higgs boson in all decay channels available 
^■s^ ■ below the top-antitop threshold is observed to be a Gaussian distribution of the Higgs boson mass with a max- 

imum centred around m u « 125 GeV, i.e. in the region of masses where a new boson has been discovered at 
f"^ , the Large Hadron Collider. This fact places the Higgs particle at a mass of "maximum opportunity' in terms of 

£N| ■ the study of its decays and couplings to other particles. Such an intriguing observation is seemingly driven by 

. p/ the different ra H -power dependence of the Higgs decay-widths into gauge bosons and fermions featuring a steep 

anticorrelated evolution below the WW decay threshold in the transition region of masses ra H a m w - 2/nw- 
TZj ■ Speculative consequences of taking this seemingly accidental feature as indicative of an (unknown) underlying 

■ dynamics of the SM are also discussed. 
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dl I. INTRODUCTION 

Oh 

In the standard model (SM) of particle physics, the masses of the electroweak W and Z vector bosons are 
i ' generated via the Higgs-Brout-Englert-Guralnik-Hagen-Kibble mechanism [Q]] which also predicts the existence of 

a new elementary neutral scalar (Higgs) boson. The fermion masses are also generated through Yukawa interactions 
with the Higgs field. The mass of the Higgs boson itself, given by the square-root of its (unknown) self-coupling 
parameter A {m^ = V2/1/Gf with Fermi constant Gf ~ 1 .16637- 10 5 GeV -2 ) is however a free parameter of the 
theory. Theoretically, the value of m u is constrained from above by requirements of unitarity of the elastic WlWl 
cross sections J2j] as well as by triviality bounds for large self-couplings A JH]. Lower limits on m K are imposed 
by vacuum stability arguments (A cannot be too small since otherwise the top-quark Yukawa coupling would 
make it negative) yfl. The recent experimental observation at the Large Hadron Collider (LHC) of a new boson at 
m H » 125 GeV in various decay modes consistent with the hypothesized SM Higgs boson |0,|5t], constitutes a ex- 
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_ perimental breakthrough in our understanding of the mechanism by which all elementary particles acquire their mass. 



Basic conservation rules imply that a neutral spinless SM Higgs boson can only decay into (i) pairs of fundamental 
fermions qq = uu, dd, ss, cc, bb, tt for quarks and t + l~ = e + e~,;U + /i~, t + t~ for leptons (neglecting neutrinos which 
are massless in the SM, see below), (ii) pairs of (virtual or real) heavy bosons VV = ZZ,WW; as well as (iii) pairs of 
massless bosons (gg,yy, as well as Zy) via heavy-particle (mostly W and top-quark) loops. The partial decays widths 
T, of the Higgs boson into the most relevant channels, i.e. the inverse decay rates to each final state kinematically 
allowed for a given value of m H , have been determined including various high-order quantum chromodynamics 
(QCD) and electroweak corrections (see g] and refs. therein). The total width is obtained by summing all the 
Higgs partial widths: r tot = 2;F/- The branching ratio of a given mode is BR = r,/r tot and, of course, the sum 
of all of them amounts to one, 2/BR, = 1. Decay widths and branching ratios to the various final-states can be 
computed using programs such as hdecay and have typical theoretical uncertainties in the range ABR « 3-10% 
for m H ss 125 GeV 131 . The SM Higgs couplings to fundamental fermions are proportional to the fermion masses nif, 
whereas its couplings to bosons are proportional to the square of their masses mL, a clear indication of the primary 
role of the Higgs particle in the breakingof electroweak symmetry. The partial ff and VV decay-widths are thus 
given by the tree-level relations (see e.g. ]8[]) 
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where N c = 3(1) for quarks (leptons), (1) 
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From these expressions one sees first that the decays to fermions are basically linear in m u but follow a 
dependence for bosons, i.e. for large enough m u the decays to W and Z are clearly preferred over fermions. [The 
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loop-mediated gg and yy,Zy decay widths - with more complex expressions not reproduced here - feature also a cu- 
bic dependence on m H , but with a smaller ri oop -prefactor than for VV decays]. Equation (O is only valid for on-shell 
vector bosons, i.e. for m H > 2my. In the region where one or both W and Z particles can be produced off-shell, the 
corresponding partial widths are smaller but rise much faster with the Higgs mass (a power-law fit to rw(m H ) in the 
range m H m mz-2mw computed with hdecay shows a dependence as fast as m 1 ^). As result, the subtreshold WW* and 
ZZ* branching ratios are larger than those for fermions already for m u as 135 GeV and 157 GeV respectively (Fig.Q]). 

In this work, we investigate whether there is a particular mass region for which the SM Higgs boson has the 
maximum number of decay modes effectively available. More precisely, we want to see which is the value (if any) 
of m H which maximizes the distribution of the product of branching ratios 

BRprodCmJ = ] ] BR H -,i(m H ) . (3) 

i 

One finds somehow surprisingly that the distribution BR pro< j(m H ) is Gaussian in the range of masses m H ss tni—lm^ 
with a maximum centred around m,. » 125 GeV and a width of cr m * 6 GeV. The reasons of such a behaviour are 
explored below. 

II. THEORETICAL SETUP 

The partial widths and branching ratios of the Higgs boson in the SM are determined following the same approach 
developed in Refs. 0H. We use hdecay (version 4.43) |@] for the fermionic, digluon, diphoton and Zy decays, and 
prophecy4f (version 2.0) fioll for the weak boson decays, hdecay includes next-to-leading order (NLO) corrections 
(or higher-order ones in some cases) for the QCD and electroweak contributions for all decay channels, whereas 
prophecy4f takes into account complete NLO and electroweak corrections as well as interferences for all final-states 
in the W,Z decays. Both calculations are then combined by modifying the Higgs total width obtained with hdecay 
according to T H = F™ - T™ - r™, + r P ™ ph 0], where r HD , T™? and T™, stand respectively for the total, and 

partial widths to ZZ and WW calculated with hdecay, and r^ oph represents the partial width of H — > 4f, including 
interferences, calculated with prophecy4f. The results of prophecy4f modify by about +10% the branching ratios 
computed with hdecay alone mostly in the mass region m H = otz-2otw- For the input parameters of the calculations 
the latest PDG values ifTHl are employed (TableUJi, with the uncertainties quoted for a few parameters being the same 
as those used in 0,Ht]. The strange-quark mass used is the MS value at two different scales Q — 1,2 GeV, and the 
charm and bottom quark masses correspond to the 1-loop pole values which are less sensitive to a s uncertainties. 



TABLE I: Input parameters used for the Higgs boson branching ratios (BR) calculations carried out in this work. 



m s (GeV) 


0.1(2 GeV),0.19(l GeV) 


m e (MeV) 


0.510998928 


mw (GeV) 


80.385 


«(0) 


1/137.035999074 


m c (GeV) 


1.28 ±0.03 


m,, (MeV) 


105.658367 


m z (GeV) 


91.1876 


G F (GeV" 2 ) 


1.1663787 • 10~ 5 


m„ (GeV) 


4.16 ±0.06 


m T (MeV) 


1776.82 


Tw (GeV) 


2.085 


Q!s(ot|) 


0.1184 ± 0.002 


m, (GeV) 


173.5 ±2.5 






Lz (GeV) 


2.4952 







The Higgs boson decays into the first-generation fermions (uu, dd, e + e~) are not directly implemented in hdecay 
since, due to their small masses compared to m H they have extremely reduced branching ratios and are of no exper- 
imental relevance. From Eq. (HJ the dependence of the first-family branching-fractions on m H can be obtained from 
the corresponding widths for heavier (e.g. second-generation) fermions via 

BRuaddC^) * {—) -BRsiK), BR e+e -(m H ) » (^-) • BR^rnJ . (4) 

The branching ratios for the quarks of the first generation are obtained using this expression with 
OT u,d,s ~ 3,5, 100 MeV. Although the resulting uu, dd (e + e~) decays are suppressed by factors of order 10~ 3 (10 -5 ) 
compared to the ss (/v + /i ) channel, the dependence of their BR on m H is the same as for other heavier fermions. 
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Since the light-quark masses are not precisely known, there are potentially large uncertainties in their associated 
BR values. As a cross-check we computed the branching ratios using the two values for the s-quark mass quoted 
in Table U and found that although the corresponding absolute BR llQd( j s§ values changed by a factor of two, the m H - 
dependence of BR ufl dd sS (m H ) is the same for any Higgs mass value. Similarly, the same changes of m UjdjS modify the 
total Higgs decay width only at the permille level given the smallness of these suppressed branching ratios. Since 
we are basically interested in the shapes of BR(m n ), and since unc ertainties in BR uds propagate just as an overall 
normalization factor of BR prod (m H ), we do not take into account the u,d,s mass uncertainties hereafter. 

III. RESULTS 

The branching ratios for all kinematically allowed channels of the SM Higgs boson are computed for masses 
m H ai 0-1 TeV using the setup discussed in the previous Section. Figure Q] shows the nine fermionic, two bosonic, 
and three loop-mediated BR in the mass range m n = 10-1000 GeV. Several observations are worth pointing out. The 
eight lightest fermionic channels (plotted here well above their respective thresholds at m K = 2m { - { ) have a relatively 
fiat BR as a function of m n up to about ~1 10 GeV, where the very fast-rising VV* decays start to play an increasing 
role which subsequently make all the ff branching fractions drop equally steeply. Beyond 2my the bosonic decays 
are roughly constant as a function of m H except for the relatively small modulation due to the onset of the tt mode 
around m K * 2m t , while all other fermionic decays decrease monotonically. The loop-mediated bosonic channels 
(gg,yy,Zy) show first a rise with m H peaking at m H as 120-140 GeV followed also by a monotonic falloff. It is 
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m H (GeV) 

FIG. 1 : Decay branching fractions of the Higgs boson into all fermion and boson pairs of the SM, as a function of the Higgs 
boson mass, computed with hdecay and prophecy4f with the parameters listed in Tableland using Eq. © for the first-generation 
fermions. The dashed vertical line indicates the position of the 125-GeV boson resonance observed at the LHC fld]. 

clear from Fig.Q]that all the "action" occurs over the mass region 2«iw) centered at around m H a: 125 GeV 
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(dashed line in Fig. [T}. The value of the W mass and the preferred decay of the SM Higgs boson into weak bosons 
are responsible for all the features seen in the distribution of branching ratios over this distinctive transition mass 
region. Indeed, the fast rise of the bosonic decays, the consequent fast drop of all (light) fermionic modes, and the 
maxima of the loop-mediated decays around the region of the newly observed LHC boson resonance, are all just 
driven by the particular value of the W boson mass 1 . 

Figure [2] shows the product of all computed branching ratios as a function of m u , Eq. 0, below the tt threshold 2 
plotted in logarithmic (left) and linear (right) scales. The band/boxes around the curves indicate the corresponding 
theoretical and parametric uncertainties, estimated as discussed in ]7[], propagating only those which are uncorrelated 
across the individual decay channels. The absolute amplitude of BR pro[ j(m H ) has no particular meaning and the 
distribution has been normalized by its integral. We note that below mz, the Zy-decay is actually null in hdecay since 
the code does not include the Z*y channel 3 . This fact explains the sudden drop seen in BR pl0C j(m H ) around 90 GeV. 
The resulting distribution can be very well fitted by a Gaussian function in the transition region from the fermion- 
to the boson-dominated regime m H = mz-2mw with goodness-of-fit per degree- of-freedom: * 2 /ndf = 61.6/68. 
Surprisingly, the distribution has not only a Gaussian shape but it has a maximum centered at m n as 125 GeV (and a 
width of about cr,,^ = 6 GeV) i.e. in the range of masses where the Higgs-like particle has been discovered at the 

LHC: m™ s = 125.3 + 0.6 H and m* TLAS = 126.0 + 0.6 0] (hatched grey band in Fig. 0. Although the fit quality 
(X 2 /ndf « 1) is very good within 3ct„, h , the product of branching ratios starts to deviate from a Gaussian in the tail 
beyond m B ~ 155 GeV. This disagreement may be corrected if one adds the tt Higgs decays including the 4-body 
contributions from diagrams involving two virtual top-quarks which are very small in amplitude (although steeply 
rising with m H ) and currently neglected in hdecay. 




FIG. 2: Product of all decay branching ratios of the SM Higgs particle into fermion and boson pairs below the tt threshold as 
a function of the Higgs boson mass, Eq. $3$ normalized to its integral. The distribution is shown in log (left) and linear (right) 
scales and fitted to a Gaussian with the parameters listed in the legends. The uncorrelated uncertainties obtained in 01 are shown 
as boxes/band. The hatched vertical band indicates the position of the observed boson resonance at the LHC UQ]. 



We have also checked what happens with BR pro d(/« H ) if one excludes one or various decay channels. As a 
matter of fact, the product of the two dominant fermion (bb) and bosonic (WW) decay branching ratios is already 



In the SM, the Z and W masses are linked through the Weinberg mixing angle. 

2 The tt branching ratio is actually set to zero in hdecay below m H » 250 GeV and is omitted in the product of BR(m H ) distributions. 

3 Likewise, hdecay does not include the Dalitz y*y mode with y* decaying into pairs of fermions - amounting to about 10% of the y y decay ratio 
according to fl2ll - which nonetheless has the same m H -dependence as BR ry and thus does not affect the shape of the obtained BRp ro[ j(m H ) 
distribution. 
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a distribution with Gaussian-like shape and a maximum at m H a 135 GeV with larger width cr,,,^ =s 14 GeV. The 
product of the three loop-mediated branching fractions alone is also normally distributed around m H as 131 GeV 
with a width cr m ~ 1 1 GeV. The consecutive addition of extra fermion-pair BR displaces the Gaussian to lower m H 
values and reduces its associated width, resulting in the final overall distribution plotted in Fig. [2] 

The first pragmatical consequence of the observation plotted in Fig.[2]is that a SM Higgs boson around 125 GeV 
has the largest possible number of decay modes concurrently open to potential experimental study. No other Higgs 
mass value has theoretically a better "signal strength" for the whole set of kinematically-allowed decay channels 
simultaneously. This fact places the Higgs boson in a mass window of "maximum opportunity" to experimentally 
study its couplings to all the gauge bosons and fermions at the LHC as well as at any other future collider. 

A priori, there is no physical reason why the product of the many (or a few) branching-ratios distributions plotted 
in Fig. Q] should be normally distributed. However, in statistics theory it is well-known that probability densities 
such as the beta distribution, P(x) oc x™(l - x)" for < x < 1, and the gamma distribution, P(x) oc x m exp(-.x) for 
x > 0, converge to Gaussian distributions 4 for large values of the powers m and n. From the tree-level Higgs-mass 
dependencies of the decay widths given by Eqs. ([D and ©: T f - f oc m iV Fyv K m\ (and Fi oop oc m^), the product 
of branching ratios for the fermionic and bosonic (and loop-mediated) decays is indeed a expression which has the 
following generic parametric form 

BR pr0 dK) oc ( P| r,<m H )) x ( J Ti(mjj " oc (mjr . m ^ + ^ x ^ . ^ . ^ + R . ^ + ^ . m 3)"« > 

i i 

for a total number of decay modes n - « f f- + «vv + »ioop = 13 (below the tt threshold) and constants A and B. Taking 
n f f = 8, «vv = 2 and «i oop = 3 one obtains a beta-type distribution with large exponents 

BRp ro dK) oc m^x(A' +B' -miy 13 (5) 

which can be indeed a Gaussian-like function for appropriate positive A' and B' constants. 

Although the observed Gaussian shape of BR pro d(f7J n ) seems to be accidental and it happens to peak around m H =s 
1 .5mw for the reasons described above, one can speculate about the possibility that there could be an (unknown) 
underlying dynamical reason that makes the product of SM Higgs boson branching fractions to peak at m n and, 
therefore, use this fact to place limits or constraints on new possible particles and/or Higgs decay channels. Making 
a simple analogy with entropy arguments in the decay of an excited system in equilibrium, one could naively argue 
that given that the physical realization of the Higgs field - which gives mass to the rest of elementary particles - is 
an unstable particle, the value of m H must be such that the Higgs boson has the largest number of ways in which its 
available mass can be potentially distributed into the known fundamental particles. Thus, if the 125-GeV resonance 
observed at the LHC is indeed the SM Higgs boson and the maximum of the product of its branching-ratios at m H 
derives from some more fundamental SM ground one could constrain the Higgs "invisible" decays into e.g. neutrinos 
or axions lfl3ll . Adding, for example, three new Higgs boson decays into (Dirac) neutrino-antineutrino pairs with 
Yukawa couplings as for the other SM fermions, i.e. via Eq. © using m v instead of m e , displaces the peak of 
the distribution to about m H « 122 GeV. It would be interesting to test the shape of the corresponding BR pl0 d(m H ) 
distribution taking into account different neutrino mass realizations e.g. in the context of seesaw mechanisms. 
Similarly, we defer for a coming study to check whether the m H -peak in the product of decay fractions present in the 
SM appears also for the (unconstrained) mass of other Hi ggs particles in the context of supersymmetric (SUSY), or 
more generally two-Higgs-doublet, extensions of the SM MM . 

IV. SUMMARY 

The fourteen branching ratios of the standard model Higgs boson into the known fermion and boson pairs have 
been computed as a function of the Higgs boson mass in the range m H « 0-1 TeV with the hdecay and prophecy4f 



4 With means fi = m/(m + n), m and widths <x = (m ■ n)/(m + n) 3 , m for the beta and gamma distributions respectively. 
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programs. The product of the branching ratios in all decay channels available below the top-antitop threshold is 
observed to be a Gaussian distribution of the Higgs boson mass with a maximum centred around m H as 125 GeV, 
i.e. in the region of masses where a new boson resonance has been discovered at the LHC. We have argued that 
such an intriguing feature is likely the accidental result of the fact that the fermionic 5 , bosonic and loop-mediated 
Higgs decay branching ratios feature all strongly anti-correlated changes around m B w 1 .5»i\y due to their different 
m u -power dependencies. As a result, the product of their branching ratios has a beta-type probability distribution 
which converges into a Gaussian for large number of decay channels. If, on the contrary, this observation is due to 
an (unknown) underlying physical mechanism - i.e. if there is any fundamental reason why the value of m H must be 
such that it has the largest number of ways in which the available mass can be potentially distributed into the known 
fundamental particles - such an effect could be used to impose extra constraints on new possible particles and/or 
Higgs decay channels in the SM and its extensions. 

In any case, the observation reported here places the Higgs boson in a mass window of "maximum opportunity" to 
experimentally study its couplings to all the gauge bosons and fermions, and therefore impose maximum constraints 
to any physics beyond the standard Model, at the LHC as well as at any other future collider. 
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